Topological transition of Dirac points in a microwave experiment 
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By means of a microwave tight-binding analogue experiment of a graphene-like lattice, we observe 
a topological transition between a phase with a point-like band gap characteristic of massless Dirac 
fermions and a gapped phase. By applying a controlled anisotropy on the structure, we investigate 
the transition directly via density of states measurements. The wave function associated with each 
eigenvalue is mapped and reveals new states at the Dirac point, localized on the armchair edges. 
We find that with increasing anisotropy, these new states are more and more localized at the edges. 
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Introduction. — Recently discovered condensed matter 
systems, such as graphene [1] or topological insulators Q, 
constitute an ideal playground to investigate the physics 
of massless Dirac fermions until now restricted to high en- 
ergy physics. Indeed, the relativistic spectrum emerges at 
the conical intersection points, the so-called Dirac points, 
in the low energy electronic dispersion relation. The large 
potential of technological applications [3] depends cru- 
cially on the properties of these quasi-particles. The ma- 
nipulation of the Dirac points by external parameters - 
from creation to annihilation - have recently raised a 
stronginterest both theoretically 0-0 an d experimen- 
tally [l, [H . By controlling the anisotropy of a honey- 
comb lattice, one can in principle move the Dirac points 
up to a transition where they merge and annihilate each 
other. This is a topological transition since both Dirac 
points are characterized by opposite topological numbers 
(opposite Berry phases) which annihilate at the transi- 
tion. If the anisotropy is increased further, a band-gap 
opens in the dispersion relation. This transition from 
a gapless (Dirac) phase to a gapped phase corresponds 
to a Lifshitz phase transition from a semi-metallic to an 
insulating phase [4-7]. The experimental investigation 
of this transition is difficult to perform with condensed 
matter systems, e. g. graphene under uniaxial strain [To| 
or quasi-two-dimensional organic conductors under pres- 
sure 11]. Artificial systems such as molecular graphene 



or ultracold atoms in optical lattices have been recently 
proposed to experimentally probe this topological transi- 
tion [g, 9]. While this effect has not been clearly observed 
with molecular graphene [q| , Ref. [H reports the observa- 
tion, using momentum-resolved interband transitions, of 
Dirac points merging with a Fermi gas in an anisotropic 
honeycomb lattice without six-fold symmetry. 

Alternatively, discrete photonic systems (e. g. photonic 
lattices or microwave cavities), with graphene-like struc- 
tures, have been employed as condensed-matter ana- 
logues 12 -15|. In microwave experiments, the tight- 



binding form of graphene's Hamiltonian [16| can be es- 
tablished by using a honeycomb lattice of evanescently 
coupled dielectric resonators [H], [Hj • Dirac points and 



signature of the linear dispersion relation have been ob- 
served I3-[iH- 

In this Letter, we take advantage of the high versatil- 
ity of our tight-binding microwave set-up to explore the 
topological phase transition. The anisotropy of the hon- 
eycomb lattice is controlled through the nearest-neighbor 
hopping parameters. We directly measure the density 
of states (DOS) and the wave function associated with 
each eigenvalue. We observe a transition between the 
Dirac phase and the gapped phase exactly for the ex- 
pected value of the anisotropy parameter defined in a 
tight-binding description. This is the first direct DOS 
measurement of this topological transition. Moreover the 
experimental set-up allows a very fine manipulation of 
the edges and a control of edge states which is not yet 
possible in condensed matter systems. We show that the 
anisotropy generates new states, at the Dirac frequency 
vd, localized along specific armchair edges. With increas- 
ing anisotropy, their extension into the bulk decreases. 

Experimental microwave graphene analogue. — For our 
experiments, we use a set of 222 identical coupled dielec- 
tric cylindrical resonators (5 mm height, 8 mm diameter 
and a refractive index of 6), hereafter called discs, placed 
in between two metallic plates. We establish a two- 
dimensional tight-binding regime where the electromag- 
netic field is mostly confined within the discs and spreads 
out evanescently. The experimental details are described 
in [3, EH- Via a movable loop-antenna, the reflected 
signal Sii is measured over a given frequency range at 
the disc positions r using a vectorial network analyzer. 
The bare frequency of an isolated disc (~ 6.65 GHz 
in the set-up considered here) corresponds to the on-site 
energy appearing in the tight-binding Hamiltonian. The 
coupling parameter t exp between two adjacent discs de- 
pends on the disc separation d and can be described, as 
expected for a two-dimensional system, by a modified 
Bessel function Kq (see Table. [Tj). Thus, by changing the 
distance between discs, one can change, in a controllable 
manner, the inter-site couplings. For a lattice composed 
of many discs, the set-up allows for a direct access to 
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FIG. 1. (Color online) Left panels: Schematics of the hexag- 
onal lattices, ai and a.2 define the primitive cell of the Bra- 
vais lattice, t (blue intersite links) and t' (red ones, along y 
axis) correspond to the nearest-neighbor couplings, d and d! 
are the corresponding distances. An anisotropy is introduced 
along the y axis: t' increases when d! decreases. Right panels: 
Calculated density of states (DOS) for an infinite system with 
first, second and third nearest-neighbor couplings given in the 
text. The red dashed line is the Dirac frequency vr> — ^o+3^2- 
The insets show the corresponding dispersion relation, (a)-(b) 
Regular honeycomb lattice, f5 — 1. The two Dirac points are 
distinguishable, (c)-(d) Compressed lattice at the topological 
transition, f3 = 1.8. The Dirac points merge, (e)-(f) f3 = 3.5, 
a band gap opens (gray zone). 
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FIG. 2. (Color online) Experimental DOS for (a) the regular 
graphene (f3 = 1), (b) and (c) anisotropic structures with f3 — 
1.8 and fi — 3.5 respectively. The dashed red lines indicate 
the Dirac frequency vr>- The gray zone shows the band gap 
with edge states at vr>- 



the local density of states (LDOS) through the measured 
quantity Sn [3, Ell- The DOS is obtained by averag- 
ing over all positions r. Moreover, our set-up allows to 
visualize the wave function associated with each eigen- 
frequency [I?}- 

Tight-binding description. — In order to interpret quan- 
titatively our experimental results, we present a tight- 
binding model on an infinite hexagonal lattice with first, 
second and third nearest-neighbor couplings. The on-site 
"energy" is given by the bare frequency The hop- 
ping amplitudes (in units of frequency) between nearest- 
neighbor sites are called t f along the ?/-axis and t along 
the other directions (respectively red and blue links on 
Fig. []Ja). The next (next) nearest-neighbor couplings 
are respectively denoted as t2 and £3. The anisotropy of 
the lattice can be simply taken into account via the first 
nearest-neighbor coupling (i.e. £'), higher-order terms 



TABLE I. Experimental anisotropy parameter ft — £e xp Aexp- 
For the regular lattice (i.e. /3 = 1), d = 15 mm and texp — 
0.016 GHz. The intersite distances d! and the corresponding 
nearest-neighbor coupling values /3t exp are also indicated. 
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15.0 


0.015 
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13.9 


0.020 


1.3 


13.0 


0.027 


1.8 


12.6 


0.030 


2 


12.0 


0.037 


2.5 


11.0 


0.053 


3.5 



a Obtained according to the Bessel fit t e xp(d) = a \ Ko('yd/2)\ 2 
+5. a = 1.9545 GHz, 7 = 0.3167 mm- 1 and 6 = 0.0049 GHz 
have been experimentally determined fl7| . 
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are assumed to be unchanged (i. e. t 2 = t' 2 and £3 = t' 3 ). 
In the basis of the two inequivalent sites forming the ele- 
mentary cell, the Hamiltonian can be written in the Bloch 
representation: 



^o + / 2 (k) /(k) + / 3 (k) 
/*(k) + / 3 *(k) ^ + / 2 (k) 



(1) 



where / (resp. / 2 and ^3) is the first (resp. second 
and third) nearest-neighbor contribution: /(k) = —(t' + 
te ik. ai +te ik-a 2 ^ respectively / 2 (k) = -2£ 2 [cosk-ai + 
cosk • a 2 +cosk • (a x - a 2 )] and / 3 (k) = -£ 3 [e* k ' (ai+a2) + 
e *k-(ai-a 2 ) ^ e *k(a 2 -ai)j frequency spectrum is given 

by z/(k) = ^0 + /2(k) =b |/(k) + /s(k)| and consists of two 
bands touching at the Dirac frequency up = v$ + 3t 2 . 
The corresponding DOS are represented in Fig. [T] for dif- 
ferent values of the anisotropy parameter f3 = t'/t. For 
/3 = 1, two Dirac cones exist in /c-space at k = = 
=b(a* — a 2 )/3 where a* are reciprocal lattice vectors. For 
larger /?, the Dirac points moves along k x axis. A topo- 
logical transition occurs at /3 C = 2 — 3^/^, where the two 
Dirac points coalesce [inset Fig.QJd)]. With the values of 
the coupling parameters extracted from the experiments 
(see below), we find j3 c = 1.8. For j3 larger than this 
critical value, the condition z/(K^) = cannot be ful- 
filled: a gap opens [inset Fig. QJf)]. As a consequence, 
the density of states is strongly affected. The typical lin- 
ear relation at the vicinity of the Dirac point (red dashed 
line) observed for f3 = 1 [Fig. QJbYl turns to a square root 
law at the transition [Fig. QJd)] 0, 0]- Then a gap opens 
for f3 > j3 c [gray zone in Fig.QJf)]. Note that next (next) 
nearest-neighbor couplings lead to (i) a lower value of (3 C 
for the transition (1.8 instead of 2 expected when only 
the first nearest-neighbor coupling is taken into account 
0,E|), (H) an asymmetry of the bands, (Hi) a slight shift 
of vb towards lower frequency and (iv) a displacement 
of the two logarithmic singularities. 

Topological transition and gap opening. — 
Experimentally, we first realize a regular graphene 
lattice with an hexagonal shape which has only armchair 
boundaries in order to eliminate edge states [see inset 
Fig. El^a)]. We gradually add an anisotropy by increasing 
the coupling t' eKp along the y axis - i.e. by decreasing 
the separation d! between the discs along y (insets in 
Fig. [2j). As detailed in the Tab. HJ the anisotropy pa- 
rameter f3 varies from 1 to 3.5. The DOS corresponding 
respectively to j3 = 1,1.8,3.5 are depicted in Fig. [21 
For the regular graphene, for which d = d! = 15 mm, 
the non-symmetric band structure and the point-like 
gap at the Dirac frequency vb are observed [Fig. [2(a)] . 
Here, due to the finite size of the system, the expected 
linear behavior of the DOS near the Dirac point appears 
only as a dip in the histogram. Moreover, the two 
logarithmic singularities have their signature in this 
finite sample (Fig. [2]). By adjusting the calculated DOS 
shown in Fig. [1(b) to the experimental one [171], the 



following values for the three tight-binding couplings 
parameters are obtained: t = 0.016 GHz, tijt = —0.091 
and ts/t = 0.071. This value of t is in agreement with 
the Bessel fit of the experimental coupling (see Tab.|T]). 
Fig. [2(b) and (c) show the evolution of the DOS when the 
lattice is compressed. For large values of /?, a gap opens 
as expected but additional states appear around vp. 
We will provide a description of these states, identified 
as edge states, in the last part of the Letter. For an in- 
spection of the topological transition, we eliminate their 
contributions by averaging the signal over bulk sites only. 
In Fig. [3(a), we plot the measured DOS obtained for six 
different values of f3 using a white-blue color scale. The 
dashed red line corresponds to the Dirac frequency and 
the gray zone is used here to improve the visualization 
of the band gap opening. These experimental results 
are successfully compared with the DOS calculated 
for the corresponding infinite systems [Fig. [3(b)]. We 
clearly distinguish the topological transition between 
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FIG. 3. (Color online) (a) Experimental DOS for /3 varying 
from 1 to 3.5. To remove edges states and highlight the band 
gap opening, the DOS have been averaged over bulk sites only. 
The arrows at f5 — 1, 1.8 and 3.5 indicate the DOS plotted 
respectively in Fig. [2(a), (b) and (c), where the average was 
performed over all sites, (b) Calculated DOS for /3 varying 
from 1 to 3.5 for an infinite system. The arrows at /3 = 1, 1.8 
and 3.5 indicate the DOS plotted in Fig. [1(b), (d) and (f) 
respectively. In (a) and (b), the gray zone shows the band 
gap opening, the dotted lines denote the linear dependence of 
the band edges and the red dashed line indicates vd- 
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FIG. 4. (Color online) Experimental wave function distributions at the Dirac frequency vd (red dashed line in Fig. [3]) for f3 
ranging from 1.3 to 3.5. 



a gapless phase and a gapped phase which occurs for 
P c = 1.8, as expected from the theoretical expression 
/3 C = 2 — 3^3 jt and the tight-binding parameters given 
above. Before the transition, the spectrum exhibits a 
point-like gap. At the transition, since the DOS close to 
the Dirac point is supposed to switch from a linear to a 
square root behavior 0, Ej, the number of states close 
to the Dirac point is expected to increase. Due to the 
finite size of the lattice, a small gap always exists, but a 
higher concentration of states at ud is clearly seen which 
reveals the topological phase transition. Finally, above 
the transition a gap opens. The gap and the bandwidths 
linearly increase with P (respectively the gray zone and 
the dotted lines in Fig. [3j). 

Edge states. — We now focus on the states which ap- 
pear at the Dirac frequency ud for anisotropic lattices 
(see Fig. [2]). The spatial distribution of these states is 
depicted in Fig. 2] where it appears clearly that they are 
indeed edge states, which may appear a priori surpris- 
ing since armchair boundaries are not expected to sup- 
port edge states [18]. Their structure will be carefully 
studied in a future work, but several remarks are of im- 
mediate interest in this Letter, (i) They appear only 
along the edges which are not parallel to the anisotropy 
axis, (ii) Their localization along the edge increases when 
P increases, (Hi) Their existence is not related to the 
topological transition observed in Fig. [3j they appear as 
soon as p > 1. These features are in agreement with 
the prediction for the existence of armchair edge states 
in deformed structures [l9[. Moreover, we find that (i) 
the intensity on one triangular sublattice stays zero, (ii) 
the intensity on the other sublattice decreases roughly as 
(1/P) 2r where r is the distance to the edge in units of 
the lattice parameter. A more extensive investigation of 
these states, as well as of the states along zig-zag and 
bearded edges in anisotropic structures is in progress. 

Conclusion. — To conclude, we have exploited the high 
flexibility of a microwave two-dimensional analogue of 
strained graphene to experimentally demonstrate the ex- 



istence of a topological transition between a gapless phase 
and a gapped phase when the anisotropy increases. We 
have investigated the spectral properties of this struc- 
ture by directly measuring the density of electromagnetic 
modes. An appropriate description of this system re- 
quires to take into account not only first but also next 
(next) nearest-neighbor couplings. The wave functions 
associated with the states appearing at the Dirac fre- 
quency are easily measurable with our set-up and they re- 
veal the presence of edge states along the armchair edges. 
They become more and more localized along the edges as 
the anisotropy is increased. Finally, the versatility of the 
set-up allows to mimic many of the spectacular properties 
of graphene and to reveal many other physical phenom- 
ena such as the influence of pseudo-magnetic field and 
associated Landau levels induced by a specific strain in 
graphene-like lattices [20| . 
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